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INTRODUCTION 


As shoMn by Gauss In 1839 the potential of the geomagnetic field can 
be represented by a spherical harmonic series of the form: 

NMAXl n n+1 „ 

V « a I i (^) [g cos m* + h sin m*] P (cos e) 
n«l m=0 

( 1 ) 

NMAX2 n ^ n „ 

+ i I i (-j) [dp cos + Sj^"* sin m^] (cos 0) 
n“l m*0 

where: a Is the mean radius of the earth, 

r, 0, ♦ are the standard spherical coordinates, and 

P^”* (In "iTOdern" methodology) are the Schmidt quasi -normal 1 zed 
form of associated Legendre functions. 

The magnetic field Is then given by: 

I , -VV . (B^, Bg. B^). (2) 

Theoretically, (1) applies exactly, at a given time, only when NMAXl 
and NMAX2 go to Infinity, under the assumption that the region under 
consideration, a < r < b , say. Is source-free. The source-free 
assumption holds nearly exactly between the earth's surface and the 
Ionosphere, but near-earth spacecraft pass through a region of "field 
aligned" currents In the auroral belt. The geometry is such that the 
field magnitude and vertical component are relaively unaffected (Langel , 
1974) but the horizontal components may have several hundred nT 
(nanotesla) contribution from these currents. This must be accounted for 
in deriving (1). In practice, the values of NMAXl and NMAX2 are limited 
by the data accuracy, by finite computer capabilities and, for NMAX2, the 
nature of Its temporal variability. The data accuracy aspect will he 
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f1+l 

discussed in a later section. In (1) the terms in (a/r) describe 
sources within r < a , or "Internal" sources, and the terms In (r/a)” 
describe sources outside r ^ b, or "external" sources. Field measure- 
ments are used to derive the coefficients h^"*, and s^*", usually 
by some form of least squares procedure. 

Both the Internal and h^*") and external and s^"’) 

coefficients are known to vary with t1n». To date, temporal variations 
in the external terms have not been included In models. Variations in 
the internal field have been modeled by expanding the coefficients In 
Taylor series In time, e.g.: 






+ •< 


(3) 


Most models include only the constant and first derivative (secular vari- 
ation) terms, although some more recent models have Incorporated the 
second derivative (secular acceleration) also (e.g., Cain et al , 1967; 
Barraclough and Malln, 1979). It should be noted that in some contexts 
the derivatives In equation (3) are combined with the factorials to 
produce the total coefficients which multiply the powers of time in the 
power series. 


The principal sources of data for main field modeling have been (1) 
permanent magnetic observatories, (2) repeat measurements at selected 
sites with Intervals between measurement of one to six years, (3) surveys 
from aircraft and ship, and (4) satellite measurements. Only the 
satellite surveys are truly global. Relevant surveys from which data are 
generally available were conducted by the Cosmos 49 spacecraft in October 
and November of 1964, by the 060-2, -4, and -6 (P060) spacecraft from 
October 1965 through July 1971, and most recently, by the MA6SAT 
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spacecraft from November 1979 tnrough June 1980. The Cosmos and P060 
satellites measured only the field magnitude, which Introduces an 
ambiguity In resulting spherical harmonic analyses (Backus, 1970; Hurwitz 
and Knapp, 1974; Stern and Bredekamp, 1975; Stern et a1, 1980). The 
permanent magnetic observatories must still be regarded as the primary 
source of Information regarding the temporal changes. Unfortunately, 
these data, and all data except the satellite data, are highly 
"contaminated" by fields originating in the crust of the earth. These 
anomaly fields can be tens to thousands of nT In any of the components 
and represent a large nc £e source when attempting to model the bulk of 
the geomagnetic field which originates in the earth's core. 

This paper describes an attempt to utilize the observatory data in a 
more optimal way by Incorporating an estimation of individual observatory 
anomaly fields Into the solution and, in so doing, to allow thl^ 
derivation of a more accurate model of temporal variation. The results 
were presented as a possible contribution to the 1980 version of the IfiRF 
(International Geomagnetic Reference Field) and the definitive 
Geomagnetic Reference Fields for 1965, 1970 and 1975 at the 1981 assembly 
of the International Association of Geomagnetism and Aeronomy (lAGA). 

MODEL DEGREE AND ORDER 

In (1) "n" Is the degree and "m" the order of any given term. 

NMAXl is the maximum degree and order for internal terms and NMAX2 for 
external terms. Gauss' original model did not include external terms and 
used an NMAXl of four because he concluded that the available data did 
not warrant the inclusion of further terms. In the years since Gauss the 
available data base has Improved considerably and, accordingly, the 
degree/order of published models has Increased. Malin and Pocock (1969) 
analyzed the question of the appropriate degree/order for models based on 
magnetic observatory data. Using data from 180 observatories they com- 
puted models from degree/order two to ten and compared the rms 
residuals. Their results are shown In Table 1. 
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TABLE 1 


Degree/Order 


rois residual (nT) 


2 

3 

4 

5 

6 

7 

8 
9 

10 


3495.7 

2218.2 

889.0 

538.7 

388.5 

364.5 

359.5 
354.3 
351.9 


They note a rapidly decreasing rms from 2nd to 6th order after which the 
decrease becomes very slow; this leveling off Is attributed to the 
"crustal noise" in the data. They conclude that "at least six orders 
should be evaluated if the core field is to be fitted within 0.5 per- 
cent". Subsequent models not using satellite data are In reasonable 
agreement with their results and most of these models were of degree/ 
order eight, nine or ten. 


The situation changed drastically with the advent of satellite data, 
particularly the surveys by the POGO satellites from 1965 to 1971. These 
data are far enough above the earth that the crustal anomaly fields are 
less than 20 nT maximum, with rms below 5 nT. Further, the coverage is 
now truly global with no large gaps. The latest published models which 
include these data are mainly of degree/order 12 (Rarraclough et al , 

1975; Peddle and Fabiano, 1976; Barker et al , 1981). Those of us who 
have been attempting to study crustal magnetic anomalies have been using 
a model of degree/order 13 (Langel et al , 1980a). This choice is 
confirmed by Langel and Estes (1982) who derived the spatial power 
spectrum (Lowes, 1966, 1974), 


■ (n*l) I 

m==o 


( 4 ) 


5 
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for n-1 to 23 using Magsat data. The derived spectrum sho»i«d a clear 
break near n>14 which was Interpreted to mean that the core field 
dominates for n<13 and the crustal field for n>15. 

The situation for the temporal derivative terms Is much more 
complicated. In the model to be presented, we have Included first 
derivatives to degree/order thirteen, second derivatives to degree/order 
six, and third derivates to degree/order four. For each coefficient we 
have calculated the ratio of coefficient magnitude to the standard error 
of the coefficient. As a rule of thumb, we assuTO that If this ratio Is 
22, the coefficient has statistical significance above the 95% confidence 
level (see, e.g. Rarraclough and Malln, 1979). 

Table 2 shows the ratio distribution for these terms, (in this basis 
we are justified In Including some terms from each of the 13 degrees of 
the first derivative and, perhaps, should have extended the second and 
third derivative terms to still higher degree/order. Investigation of 
this question Is beyond the scope of the present paper. 

Aside from the question of the accuracy to which Individual 
coefficients are determined, one has to ask what the descriptive and 
predictive properties are for models with higher derivatives. We will 
address this question In a subsequent section. 

There remains the question of external fields (NMAX2). As Malin and 
Pocock (1969 and references therein) point out, various attempts have 
been made to separate the Internal and external parts of the field. 
However, the results have been widely different; in some cases 
unrealistically large external fields have been found, but have not been 
statistically significant. Again, a change has come with the global 
vector survey by MAGSAT. Langel et al (1980b) used data from November 
5-6, 1979 to derive the MGST(6/80) model. This model Included 
statistically significant external terms of degree/order equal to one. 
Furthermore, plots of the MAGSAT data clearly demonstrated the need for 
such terms to adequately represent the data. However, it Is known that 
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mu 2 : 

Sumnary of Distribution of Ratio of Coafficent 
Magnitude to Coafficent Standard Error. 


Nufl^ar 

Daorae of Terms 


RATIOS 

>100 10-100 2-10 <2 


1 

3 

2 

5 

3 

7 

4 

9 

5 

11 

6 

13 

7 

15 

8 

17 

9 

19 

10 

21 

11 

23 

12 

25 

13 

27 


First Derivatives 
1 
3 
0 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 


2 0 0 

0 2 0 

6 1 0 

6 3 0 

a 3 0 

10 3 0 

9 4 1 

12 4 1 

12 6 1 

8 8 5 

6 15 2 

2 17 6 

4 11 12 


1 

2 

3 

4 

5 

6 


3 

5 

7 

9 

11 

13 


Second Oerivatives 

0 

n 

n 

n 

0 

0 


3 

1 

6 

5 

7 

12 


n 

4 

1 

3 

3 

0 


0 

0 

n 

1 

1 

1 


1 

2 

3 

4 


Third Oorlvatlves 

3 0 3 

5 1 3 

7 0 7 

9 0 5 


0 0 

1 0 

n n 

3 1 
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the sources of the external fields, which are the magnetospheric ring 
current, magnetopause current and magnetotall current. var> widely with 
time In a fashion not yet amenable to this type of model and, moreover, 
vary more strongly as functions of local time than of longitude as set 
forth In equation (1). While external terms may be valuable as an 
Indication of the average quiet level of external fields. have no 
compelling reason to Include them In models describing extended periods 
of time. 

METHOD OF ANALYSIS 


The iTOthod of determining the parameters of (I) Is essentially the 
same as that described by Cain et al (1967) with revisions to Include 
higher temporal derivatives and to incorporate magnetic observatory data 
In a new way. Because of the data types Involved the problem Is 
non-linear and must be solved Iteratively. The Bayesian least squares 
estimation equations are as follows: 

‘»n*i ■ (5) 


where 

A Is the partial derivative matrix of the measure- 

ments with respect to the parameters 

p Is the vector of adjusted parameters 

6y Is the vector of residuals, I.e. measured data 

minus predicted value from previous Iteration 

W Is the weight matrix for the measurements 

flg is the a priori parameter covariance matrix 
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PQ Is the a pHori estimate of the parameters 

and the estimate at the (n<«‘i)st Iteratiwi Is 

Pn+1 ■ Pn + «Pn+l • (6) 

In the notation of Cain et a1 (1967), each measured quantity, say 
C, Is a function of the coefficients h^*”, and s^*” from equa- 

tion (1), denoted collectively by "p", and of the standard r,e,^,t 
space-time coordinates: 

C » C(p;r,e,^,t) . (7) 

The partial derivatives In A are then the set 3C/ap. The weight matrix’ W 
Is diagonal and formed from the standard deviations of the various 
measurements, so that for the 1th measurement with sl^na 

“i1 ““F • (8) 

As already noted, the data from magnetic observatories represent 
the most useful data set for determining the temporal variation of the 
internal field. The Incorporation of such data directly into a main 
Internal core field model, however, suffers from the fact that the 
magnetic field measured at the observatory may have a significant 
contribution due to local crustal fields. The field at an observatory Is 
represented as the vector sum 


I . *5l (9) 

where Is the Internal core field contribution from the scalar 
potential of equation (1) and Is the local anomaly field, which may 
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change appreciably over a distance of a few kllcvneters. The time 
derivative of Is assumed to be negligible so that 


• • 

I ■ . (10) 

A global satellite data set, on the other hand. Is comparatively 
free from the effects of crustal anomalies and Is certainly free from the 
effects of the higher-amplitude, more localized anwnalles. Assuming that 
for a main field model the crustal Influence on satellite data may he 
treated as random noise, these data may then be used In conjunction with 
observatory data to Isolate the non-core fields at the observatory. In 
terms of the algorithm of equations (5), (7) and (8) this Is accomplished 
by writing, for c^ch measured component at each observatory 

C “ C(p,C^ ;r,0 ,t) , (12) 

where Is the anomaly in that component at that observatory. Terms 
such as 9C/3C, are then added to the matrix, A , of partial derivatives. 

a 

The vector, p, of adjusted parameters from equation (5) then Includes p 
and all anwnaly components. This procedure allows the data to properly 
distribute their Influence among the temporal and constant parameters In 
a least square sense. Such a solution Is v^ll determined only when the 
satellite data, not strongly sensitive to the crustal fields, are 
Included In the solution. A forward elimination technique Is used for 
the anomaly components In accumulating the normal matrix A^WA, so the 
cost of obtaining the solution for the model coefficients Is greatly 
reduced. The values of the anomaly components, If desired, are then 
obtained by back substitution. 

DATA SET 

The data utilized In the present analysis Include data from (1) 
magnetic observatories, ( 2 ) repeat stations, (3) marine surveys, (4) the 
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PC^O satellite surveys and (5) the NAfiSAT satellite survey. The POfiO 
data consisted of tN 47(^ scalar observations used for the P(K^(8/71) 
model (Langel, 1974) augmented with 24(KK1 quiet 080-8 (Htservatlons from 
mld-1969 through early 1971. The HMISAT data set, consisting of scalar 
and vector measurements from Noven4>er 5 and 6, 1979, 1$ Identical to that 
used for the i^ST(6/80) model (Langel et a1, 1980b). The nviasurement 
standard deviation used to weight the P060 data was 7 nT and for MA8SAT 
data was 10 nT, based on fits to these data alone. 

Annual means data were taken from 148 observatories, selected on the 
basis of geographical distribution, longevity of measurement 
availability, and data quality. This '*esulted In our not using some high 
quality observatories regions of higher observatory density, moreover, 
only those annual means with three vector components we a accepted. The 
observatories utilized are listed In Table 3 together with the anomaly 
vector bias from the solution, the time of oata availability and the 
o. used in equation (8). The spellings used are those of the W)AA Itor Id- 
Wide Magnetic Survey of the National Geophysical and Solar-Terrestrial 
Data Center, In cases where a significant change at the observatory 
(e.g., shift of location, change of instrument) was known to have 
occurred, the data from that observatory were broken Into subsets which 
were treated Independently with respect to the anomaly bias vect(»rs. In 
all cases the measured elements were converted to X, V and Z components 
where X is north, V Is east and Z Is down In a geodetic coordinate system 
assuming an equatorial radius of 6378.165 km and a reciprocal flattening 
of 298.25. To determine the in Table 3 a least squares quadratic 
curve was fit to each component at each observatory and the standard 
deviation of the data to the solution adopted as o^. 

To fill in surface areas void of data, selected marine survey and 
repeat station data were utilized. In order to accommodate the non- 
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Table 3 

MEASUREMENT 

ANOMALY BIAS TIME SI&MAS 


STATION 

UT 

LONG 

T" 

“T 


INTERVAL 

x“ 

“T 

“ z 

«AOIl 

«.0S 

M.7R 

AOO.S 

IS.l 

lOO.R 

19R0.S-70.S 

7.1 

l.R 

1.0 

ACINCOIMI 

A1.M 

-70.77 

-AS.7 

17S.7 

-loe.o 

lORO.S-RO.l 

0.7 

1.0 

A.S 

M4«T 

M.SA 

-r7.ro 

S.7 

lO.R 

-lOi.7 

10R1.0-70.S 

17.7 

IS.O 

7o*0 

AbUAli 

U.RA 

77.07 

-00.7 

ORl.7 

RRO.A 

1RR0.S-7I.S 

1R.1 

7.S 

0.0 

AkM 

AA.N 

7R.07 

17S.0 

a7.7 

-1S7.1 

10r1.S-7r.S 

S.R 

10. 0 

lR.7 

M.M«U 

!*.•» 

•2e«« 

-R.O 

0.R 

77.7 

i;i-o.s-7r.s 

R.7 

• •2 

A.l 


•Ai.St 

177.77 

-AA.O 

-IR.O 

70.1 

10rC.->7R.S 

S.O 


• •ft 

AMMMkA|N*6A^ 

U.lS 

70 .RO 

710.S 

-OO.A 

-170. R 

lMO.S-70.1 

11.1 

o.s 

17.1 

*>1« 

-iA.Rt>47l.7« 

-RR.O 

200.1-1017.7 

10R0.S-7R.S 

7.R 


1.7 

AiNlU* 

M.SR 

11.17 

-0.A 

07 .7 

0.9 

10rO.S-7r.S 

0.2 


S.O 

MMNUNf liLm 

•M.RA 

-RA.7R 

RO.l 

-Ol.S 

S17.0 

19R0.S-R7.S 

l.A 

1.1 

7.S 

•MM* kAR( 

•AaSi 

-00.01 

1A0.0 

-A 1.0 

-101.0 

ioro.s-ro.s 

2.0 

•«f 

lR.7 


•A'.il 

-0R.01 

7 10. A 

-10s. 7 

-llA.l 

10R0.S-7R.S 

0..0 

10.1 

1.1 

•«MUI 

A.AA 

lO.RR 

-OS.R 

a;.o 

os.e 

IRO.S-RO.S 

f 

17.7 

S.2 



lO.SR 

-110.1 


1«7.1 

10RO.S-7A.S 

7.7 

27.0 

l.S 

•m«or 

n.s»-i>*.74 

71.0 


-17.0 

1OR0.S-V7.S 

ll.S 

0.7 

2ft«ft 

»USll 

•i.RA 

70.70 

OR.S 

1 <j.S 

llA.O 

10RO.S-7R.S 

A. 7 

0.0 

1.1 


AR.R7 

71.00 

-A71.0 

10.7 

110.0 

lORS.S-iR.S 


?.R 

11.0 

IlNi* 

-A. 77 

IS. 17 

-lAA.A 

•110.0 

-177. S 

1OR0.9-RR.S 

R.1 

^•2 

A. 4 

•JORNOT* 

7a. SO 

10.70 

-107.R 

AO. 7 

37.0 

lOOO.S-RA.S 

0.7 

l.R 

S.O 


7A.M 

10.70 

-110.0 

'/O.R 

7.S 

IOrS.S-RO.S 

0.0 

l.S 

7.1 


7a. SO 

10.70 

-177.1 

A7.1 

20.0 

1070.S-72.S 

0.7 

l.R 

9.1 

MUWMR 

aO.1A-IOS.7A 

-S7.S 

Sl.S 

-IrO.1 

10RA.S-77.S 

7.0 

1.1 

0.0 

•r«o 

-M.07-11<I.S< 

-7A.7 

17.S 

-1A7.7 

lOoi.S-RO.l 

^•9 

7.0 

77R 

CASTLi itOCA 

S7.7A-1M.1S 

-170.R 

-IR.S 

S.7 

1070.S-70.S 

0.0 

4« 0 

1.0 

CM A* 

77. SS 

101.01 

-IIO.R 

-00. S 

-2.0 

10RO.S-7S.S 

n.R 

7. A 

79.1 

CMMOM ROaiT 

AA.07 

7.7R 

-70.7 

-2A.S 

llA.O 

10O0.S-7S.S 


S.O 


CMkTUSRIli 

77.77 

10A.70 


-os.o 

*#0«i 

10RO.S-7R.S 

11.1 


21. R 

COINMA 

AO. 77 

-R.A7 

-S.A 

-72.0 

AO.O 

10rO.S-7r.S 

0.0 

7.S 

••ft 

COU.iiik 

•A.RR-1A7.AA 

-17.7 

-rO.7 

-107.1 

1000.S-77.S 

5*9 

• •3 

l.A 

0A1.AAS 

S7.«0 

-0R.7S 

-IS.l 

21.0 

-71. S 

10RA.R-7A.S 

lO.A 

7.0 

1.7 

9US0N 

7S.SS 

OO.S7 

-llA.S 

-lAO.O 

-270.1 

10RO.S-R7.S 

9*9 

S.2 

7. A 


71. »S 

00. >7 

-lOA.S 

-1A7.R 

-2A0.R 

10r1.S-7r.S 

9*0 

S.7 

?.A 

OOllAAS 

•7.07 

A. 17 

-01.7 

-70.0 

-71A.0 

10R0.S-7A.S 

S.i 

• •2 

7.7 

DOUAMS 

ra.ie 


-7.7 

-1R.0 

0«. A 

10RO.S-7R.S 


i.3 

7. A 

OWMMT tiua»U>i 


lAO.Ol 

-177.0 

-100.7-2R07.R 

10O0.S-77.S 


IS. 7 

It.l 



lAO.Ol 

-107.0 

-A10.1-28S2.A 

1071.S-7S.S 


19.7 

IR.I 

0W»MfTl 

A7.0A 

aa.70 

-710.A 


-02.7 

10RO.S-7R.S 

7.9 

3«* 

17.0 

OTNil 

SO. 77 

M.IO 

-7R.R 

00.1 

117.0 

IOra.S- 7R.S 


9.7 

1.7 

AiMT CMwaCMlLC 

SR. 77 

-OA.IO 

-1A0.7 


-7SS.0 

10rS.S-7r.S 

S.A 

• •ft 

0.0 

aatiiiatcaiAMC 

10.71 

-77.17 

20.1 

•^0.7 

120. S 

1009.S-77.S 


l.A 


aWtriMM 

S.A7 

-71.7a 


•aO.7 

21S.0 

10R0.S-7A.S 

2.0 

^•0 

71.0 

AiMSTMAikOaauca 

AO. 17 

11.70 

-11.7 

-S.A 

lA.O 

10RO.S-7R.S 



A.S 

OflAflCi>P A 

-11.70 

iis.os 

-Sl.O 

-127.0 

122. S 

10R0.S-7R.S 

2.; 

S.A 

17.7 

COiNIAtM 

00. 7 A 

-S1.S7 

707.R 

-217.1 

arO.O 

1000.S-71.J 

s.7 

l.S 

O.S 

CMMTAfliMATA 

A1.rO 

117.17 

-11.1 

-2 l.S 

-17.0 

10RO.A-7R.S 

1A.7 

7.7 

la. 2 

caAAT AMAc* a 

ii.i1 

-77.70 

7S0.1 

lll.S 

-SO.O 

100R.S-7S.S 

10.0 

««3 

IA.2 

caucAA 

AA.rI 

70.77 


-At. A 

-91.1 

10sO.S-7r.S 

S.A 

7.1 

R.O 

CAAM 

11 . so 

1AA.R7 

17A.0 

lOA.l 

-7.7 

1000.S-77.S 

o.i 


S.2 

HAALiT lAT 

-7S.S7 

-7r.ro 

-7R.1 

A00.7 

-7.0 

lORO.S-RR.S 


11.1 

IS.l 

HARTLAM 

SI. 00 

-A. AO 


7.0 


10R0.S-7R.S 

s.7 

• •• 

S.O 


>A.r1 

10.01 

11.0 

-IS 7. 2 

-71. R 

IOrO.S-rO.S 

l.J 

2.1 

9.2 

hcraanua 

-1A.A7 

10.77 

-17.R 

2A.0 

-S.7 

10R0.S-7S.S 

0.0 

2.0 

1.1 

HOMliLUkU 

71.17-lSO.OO 

-17A.A 

117.0 

-MA.R 

10O1.S-77.S 


7.7 

A.l 

MUANCAV0 

-17. OR 

-7S.1A 

00.1 

AO.S 

2AR.0 

10R0.S-7R.S 

11.0 

1. 1 

IR.O 

NTt#4tA«A0 

17.A1 

70.SS 
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observatory surface data In the solution In a consistent manner » tech- 
niques »#ere used to remove the major parts of the crustal contributions. 
From the available scalar marine data over the years 1<)70-1974, 39 long, 
straight tracks of length greater than 120(1 km Mere selected. A Iom pass 
filter Mas applied to each track, removing anomaly Mavelengths shorter 
than 500 km. Approximately 300 masuremnts Mere then taken along the 
filtered tracks. A measurement standard deviation of 10 nT Mas used In 
Melghting the data In the solution, although analysis of crossing points 
for the 39 tracks indicated differences on the order of 50 nT Mhen 
corrected for secular variation. The higher Meighting reflects our 
regard for the relative Importance of these data In an area devoid of 
other surface observations. Approximately 600 measurements from 150 
repeat stations Mere utilized to fill sparse data regions In Central and 
South An^rlca, Africa and Australia. Only stations Mith three or more 
occupations and good data quality Mere accepted. As there Mere generally 
Insufficient data available for the repeat stations to solve for 
Independent anomaly bias vectors, quadratic polynomial fits to these data 
Mere time differentiated to remove the assumed constant crustal 
Influence. Measurement standard deviations for the "differentiated" 
repeat data Mere arrived at by utilizing the "differentiated" value, 
together Mith observatory and marine data, to derive a degree/order eight 
model. The standard deviations to that model Mere 0.2 deg/year for p 
and 1 and 5 nT/year for ft, 1 and 6 , taken collectively. These Mere 
adopted as measurement standard deviations for the repeat station data In 
the present analysis. 

RESULTS 


In order to test the usefulness of solving for anomalies at the 
observatories and of Including second and third time derivatives, a model 
not Incorporating MAGSAT data Mas developed and its prediction capability 
tested by comparing It to the MRST(6/80) model based only on MAGSAT 


data. The test model, denoted PMAG(7/80), as well as the final model, 
denoted 6SFC(9/80), was of degree/order 13 In Its constant and first 
derivative terms, of degree/order six in its second derivative terms and 
of degree/order four in its third derivative terms. Table 4a gives the 
coefficient values and their first derivatives (secular variation), and 
Table 4b the second and third derivatives, for the fiSFC(9/8n) model. The 
standard error is given in parenthesis beside each coefficient. This is 
to be used as an error estimate only with caution because its accuracy 
depends upon (1) the accuracy of the used to weight the data (and the 
validity of the assumptions that the data are uncorrelated and have 
gaussian error) and (2) the validity (accuracy and completeness) of our 
model. In practice the standard errors tend to underestimate the actual 
error but are nonetheless useful as an indication of the magnitude of 
error in each coefficient and of the relative accuracies between 
coefficients. 

Table 5 summarizes a statistical evaluation of these models and two 
other recent (pre-MAGSAT) models. Model AWC75 (secular variation 
part) was derived by Peddie and Fabiano (1976) using data from 1967 
through 1974 and model WC80 was derived by Barker et al (1981) using data 
from 1950 through 1980, but their secular variation model is thought to 
be applicable mainly from 1974-1977. AWC75 and WC80 Include constant 
terms through degree/order 12 and first time derivatives through 
degree/order eight. In Table 5 MRST(6/80) is included as a standard of 
comparison for the 1980 epoch. Of the pre-MAGSAT models it is seen that 
PMAG(7/80) is the best predictive model. This indicates that, at least 
over a three-year interval, its temporal derivatives are not wildly 
varying beyond the data span. It should be noted that a similar model 
which did not include the observatory anomaly solution performed badly 
when used as a predictor, i.e, the presence of the anomaly solution 
affects the sc:>.'tion for the temporal terms in a positive way. The 
GSFC(9/80) mocjl represents the 1980 field well because it incorporates 
MAGSAT data. It is, however, slightly deteriorated from MGST(6/80). This 
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Table 4b 


Q8FC (9/BO) MAf»4ETIC FIELD MODEL 
(MEAN RADIUS OF THE EARTH IS S3712 KM; MEAN EPOCH IS 1980.01 
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rnldMls to 6SFC(9/tO) Mtrt taktn Including tha WST((i/Kn) external tamt with (ifi;C(d/gn) 
becauM the i^gsat data used to terlve (<!1FC(9/8fl) wet corrected for that external field. 

**9MM(7/n) and 6SFC(9/II0) Include ettlnatei of observatory *anorial1et‘* which ire taken into account 
when coagMitlng mlduals>^ 

—*0 o e s not Ineli^ data at latitudes greater than M* or test than >50*. 


Is not unexpected since It applies to a 20-year tlnw period »«hereas 
fiWiST(6/80) Is a direct fit to the two days of data Involved. 

Canparlson with observatories shows typically high residuals for the 
models which do not estimate observatory anomalies (AWC75 and WCRO). On 
the other hand, the observatory residuals for PMA"(7/80) and RSFC(9/80) 
are of the range one would expect from the published accuracy estimates 
of the observations, the Internal consistency of the data from Individual 
observatories, and the characteristics of untm>deled temporal variations. 
Note that no measurements with extreme residuals have been eliminated 
from this calculation, as Is often done. This is because such extreme 
values are often due to the anomaly values solved for In PMAR(7/8n) and 
GSFC(9/80) and so to compare the models all measurements should be 
utilized. 

For further comparison, and to get a quantitative measure of the 
predictive capability of the models, statistics were computed of each 
model versus observatory data on a year-by-year basis. In doing so, we 
adopted the statistical measure used by Mead (1979) In a similar 
analysis, namely, half the width of the median 6B% of the residual 
values. This is designated oT , and would be equivalent to one standard 
error, a , if the distribution were normal. The reason for adopting this 
measure is to avoid the situation where a few very large residuals 
dominate the statistics. Note that this is a different statistic than 
used in Table 5, and so the numbers are not directly comparable. In the 
plots to follow, some of the year-to-year changes in o are due to a 
changing distribution of magnetic observatories. Figures 1-3 show the 
variation of 7 with time for the X, Y, and Z components of observatory 
data relative to five field models: WC80, IRS75, AWC75, PMAR(7/8n), and 

GSFC(9/80). IGS75 is the designation for the model derived by 
Parraclough et al (1975) from all available data from 1955 through about 
1974. ^ is rounded to the nearest 5 nT. The observatorii's used in this 
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evaluation are those listed In Table 3 (I.e. those used In cterlvlng 
GSFC(9/80) and those listed In the Appendix. The statistic for 
GSFC(9/80) Is computed without taking Into account the local anomaly 
solution. 

Examination of these plots shorn that over Its "lifetime** the 
GSFC(9/80) model describes the observatory data as well as or better than 
the other models. The deterioration of Its representation prior to 1960 
Is also apparent. Of the models plotted. IGS75 has the longest useful 
"lifetime", presumably because the span of the data used to derive the 
model extended to the 1950*s. 

It Is very apparent that these models, considered collectively, 
suffer from a' data limitation such that a spherical harmonic analysis of 
reasonable degree/order cannot represent the data with a 7 of better 
than about 100 nT In X. 90 nT In Y and 150 nT In Z. We attribute this to 
the presence of "crustal noise" In the data and believe that it Is 
the fundamental limiting factor both on the accuracy of models based on 
such data and on our ability to evaluate model accuracy using such data. 

A clearer picture of the model degradation Is found In Figure 4. 

Here the statistics for GSFC(9/80) versus three components of observatory 
data are plotted taking Into account the local anomaly solutions. The "a 
for each component Is now In the 5 - 20 nT range from 1962. B on. These 
values are commensurate with the accuracy of measurement at 
observatories; I.e.. we believe we have very nearly eliminated the effect 
of crustal anomalies and that the statistic Is now dominated by the 
actual measurement noise. Deterioration of the model begins at about 
1961.5. within the da*'d Interval used by the model. This Is attributed 
to the lack of satellite data, with their global coverage, prior to 
1965. The amount of deterioration then Increases rapidly outside the 
data Interval. It Is roughly doubled at 1958.5. Comparison of Figure 4 
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with Figures 1 - 3 Indicates that the flKxtel deterioration becomes 
con^arable to the "crustal noise" at about 19S6. For mahy purposes this 
might be taken to be the useful limit of extrapolation of the (mwtelt 

Figures 5-8 show the yearly averages at a series of observatories 
together with the field predicted by 6SFC(9/80). Examination of these 
plots shows the need for the third time derivative, particularly for the 
X component at All bag, Boulder, Gornotayezhnaya and Guam, and the Y and Z 
components at Gornotayezhnaya and Guam. 

Figures 9-10 show data from two observatories utilized In 
obtaining the solution. We have not, then, calculated the crustal 
component for these stations and none is included In the plots. This Is 
most evident at Kodalkanal. Examination of Figures 9-10 shows that 
although the magnitude of the model differs from the data, the temporal 
change of the data Is well represented throughout the I960 - 1980 tlnw 
period. Furthermore, the third t1tf» derivative 1s Important for the Y 
components at Hurbanovo. 
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DISCUSSION AND CONCLUSIONS 


Observatory and other surface data are Inadequate In space and time 
for determining accurate field models. Hints of the limitations of these 
data for defining the main field have repeatedly surfaced, as for example 
In the study already cited by Malln and Pocock (1969). Comparing Figure 
1-3 with Figure 4 brings this limitation Into clearer focus. The fact Is 
that there Is a "noise" with an "rms" (o) of about 100-150 nT In the 
observatory data set. fcte attribute this to the existence of crustal 
anomalies, although It Is conceivable that error sources such as poor 
knowledge of data location. Instrument Inaccuracy or local magnetic 
contamination are contributing factors at some locations. Crustal 
anomalies are worldwide In distribution, have amplitude as high as 
several thousands of nT, and have a broad spectrum of spatial wave- 
lengths. Because of the poor spatial distribution of surface data, any 
attempt to even partially model the anomaly field with such data will be 
plagued with aliasing problems. We conclude that these data alone 
provide sufficient accuracy only for models up to a degree and order of 
about eight. Satellite data are not Immune to these limitations, 
although they are not as severe. Lange! and Estes (19fil) showed that 
crustal anomaly fields dwninate the Magsat data for wavelengths shorter 
than those corresponding to degree and order fifteen. This means that 
with present methods of modeling, we cannot determine the main field 
representation beyond degree and order thirteen or fourteen, regardless 
of the quality of our data. 

What must be appreciated Is that the limitations of the data are 
limitations not only on modeling but also on the evaluation of models. 

For example. In considering Figure 1, the differences between models A, 

I, P and G between 1960 and 1975 are at the "noise" (anomaly) level of 
the observatory data and are not likely to be truly significant. To 
generalize this conclusion, a model derived from satellite data with 
accuracies of, say, 10-20 nT, cannot be accurately evaluated using 
observatory data unless the model Is In error by several hundred nT. 
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Thi method of f<M* observatory biases or anomalies presented 

hunt offers a partial aolutiM to this (iita llnritatlon. Further work Is 
required to determine the depen(P>nce of the bias determination on other 
iWKl^l parameters such as the time span Involved, the degree and order of 
the constant and temporal terms, etc. (^e step In this direction has 
been taken in that we have determined the biases for several models of 
differing degree and order with only small changes In the values of the 
biases so determined. 

Further such tests need to be performed. Moreover, It needs to be 
determined If the calculated biases are In reasonable accord with, say, 
aeromagnetic anomaly data. This Is complicated by the fact that the 
untangling of the aeromagnetic anomalies from the main field model used 
to reduce the anomaly data Is not straightforward. 

Adequate representation and prediction of the temporal change of the 
earth's main field has. In our minds, not been achieved. In this paper, 
we have taken two steps which we believe Improve the situation. First , 
we have Included solution for the localized observatory anomaly fields 
and, second , we have utilized third time derivatives. The success of the 
second step depends upon Inclusion of the first. 

The use of second and third derivatives has certainly resulted In a 
model of good accuracy for a twenty year period. The advantages of this 
over the use of several Individual models for a shorter period of time 
are continuity and that the highly accurate satellite data are allowed to 
provide some constraint at times remote from the data epoch and so, 
hopefully. Increase the overall accuracy. Tests Indicated that 
PMAG(7/80) predicted two to three years forward In time more accurately 
than other Pre-Magsat models. Comparison with observatory data (Figures 
1-3, 5-10) Indicates that GSFC(9/80) "predicts" (backward In time from 
1960) within the 100-200 nT level for about four years. However, for 
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longer prediction, such models are clearly not suitable because the 
higher derivative terms begin to d(mi1nate and the resulting error 
Increases at a more rapid rate than for models with, say, only 
first derivatives. A solution might be to taper the higher temporal 
derivatives to zero a few years beyond the existing data interval. This, 
however. Is quite ad hoc and has little real justification. Another 
approach would be to continue to add temporal terms, including higher 
order derivatives, until, hopefully, some sort of convergence is 
achieved. From the discussion of Table 2, it is likely that higher 
degree and order second and third derivatives are significant. It is not 
clear, however, that temporal convergence will ever be achieved within 
any reasonable computer limitations or even that it is possible. The 
problem, of course, is that model constraints based on the physics of the 
core dynamo are not built into these models. In our view, until this is 
done it will not be possible to derive truly adequate "forecast" models. 
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Appendix: Observatories used for the statistics of Figures 1-3 but not 
Included In the model solution are: Abinger, Abisko, Acacias, Aso, 
Averroes, Baguio, Bouzareah, Budakeszl, Budkov, Ca^rldge Bay, 
Castellacclo, Centro fieofisico, Cheltenham, Davao, Dehra Dun, Druzhnaya, 
Ebro, Eights, El Ablod Sidl, Ellsabethvllle, Eskdalemuir, Blbllmanna, 
Gonzalez Videla, Hallett Station, Heard Island, Heiss Island, Helwan, 
Hollandia, Hurbanovo, Ibadan, Isla da Pascua, Jassy, Jullanehaab II, 
Kanoya, Kanozan, Karavia, Katuura, Kiruna, Kodalkanal, Kuyper, L. America 
III, L. America V, Lazarevo, Luanda Golf, Manhay, Maputo, Marlon Island, 
Mizusawa, Monte Capellino, Murmansk, Nagycenk, Nantes, Nitzanim, Norway 
Station, Novo-kazallnsk, Orcadas del Sur, Patrick, Pendell, Port -Alfred, 
Regensberg, Rol Baudouln, San Fernando, Simferopol, Srednikan, 

Stonyhurst, Surlarl, Swider, Taipei, Tangerang, Thule I, Tikhaya Bay, 
Toledo, Tulsa, Voroshilov, Vykhodnoy, Watheroo, Wien Auhof, Wien Kobenzi, 
Witteveen, Yellowknife. 
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Figurt 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Figure 6: 


Figure 7: 


Variation of with time for obstrvat^ data relative to 
five field models, not Including cdiservatory anomalies. 

Variation of Oy with time for observatory data relative to 
five field models, not Including observatory anomalies. 

Variation of with time for ^servator^ data relative to 
five field models, not Including <MMervatofy anomalies. 


Variation of o with time for 6SFC(9^), Including observa> 
tory anomalies. 

Co^arlsmi of Allbag annual means (X,Y,Z) to values computed 
from 6SFC(9/80), Including the observatory anomaly. Data 
from this observatory was used In the solution for GSFr(9/B0). 

Coa^)ir1son of Boulder annual means (X,Y,Z) to values cor^)uted 
from 6SFC(9/80), Including the observatory anomaly. Data 
from this ob'iervatory was used In the solution for RSFC(9/80). 

Comparison of Gomotaye^hnaya annual means (X.Y,Z) to values 
computed from GSFC(9/80), Including the observatory anomaly. 
Data from this observatory was used In the solution for 
GSFC(9/80). 
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Figurt 8: Comparison of &iam i'nnual moans (X.Y.Z) to valuos conputod 
from 6SFC(9/M), Including tlw obsorvatory anomaly, nata 
from this observatory was used In the solution for GSFC(9/80) 

Figure 9: Comparison of Hurbanovo annual means (X.Y.Z) to values 

computed from 6SFC(9/80). No observatory anwsaly was com- 
puted because data from this (^servatory ms not used In the 
solution for 6SFC(9/80). 

Figure 10: ComparISM of Kodalkanal annual means (X.Y.Z) to values 
confuted from GSFC(9/80). No observatory anomaly was 
computed because data from this observatory was not used in 
the solution for GSFC(9/80). 
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